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o: 

^vj The vortex states in a thin mesoscopic disk are investigated within the phenomenological 

Ginzburg-Landau theory in the presence of different "model" magnetic field profiles with zero av- 
erage field which may result from a ferromagnetic disk or circulating currents in a loop near the 
) ^ superconductor. We calculated the dependences of both the ground and metastable states on the 

^H . magnitude and shape of the magnetic field profile for different values of the order parameter angular 

moment, i.e. the vorticity. The regions of existence of the multi-vortex state and the giant vortex 
state are found. We analysed the phase transitions between these states and studied the contribu- 
tion from ring-shaped vortices. A new transition between different multi-vortex configurations as 
the ground state is found. Furthermore, we found a vortex state consisting of a central giant vortex 
surrounded by a collection of anti-vortices which are located in a ring around this giant vortex. The 
limit to a disk with an infinite radius, i.e. a film, will also be discussed. We also extended our results 
to "real" magnetic field profiles and to the case in which an external homogeneous magnetic field is 
present. 
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^ , I. INTRODUCTION 



Recent progress in microfabrication and measurement techniques makes it possible to study the properties of small 
superconducting structures, so called mesoscopic samples, with sizes comparable to the penetration deplli A and 
the coherence length ^. Mesoscopic disks have been one of the most popular and exciting study objectsEI E3 in this 
respect. The behavior of such mesoscopic samples in an external magnetic field is strongly influenced by the boundary 
conditions, sample size and geometry and may lead to various superconducting states and phase transitions between 

■^ \ Motivated by recent experimentsoliil, we study the properties of a superconducting disk in the presence of a step-like 

I~~«- ■ external magnetic field. The step-like field profile is a model magnetic field profile for a ferromagnetic dot or current 
^^ loop placed on top of the superconductor. These profiles have the important property that the average magnetic field 
'""' is zero. We investigate the influence of step height, step profile, and ratio between step width and radius of the disk, 

^^ on the superconducting phase diagram. 

Previous investigations of structures with naagpetic dots were limited to experiments with superconducting films 

deposited on regular arraizSjjif magnetic dotscJiiB and theoretical studies of single magnetic dots embedded in a 

thin superconducting filmllZlilj. The common problem was that, for magnetic dots, the strong field present inside 

i-A . the ferromagnet suppresses the superconducting order parameter, and in such situations it is appropriate to adopt 

^ ' a boundary condition in which the order parameter itself vanishes. This spoils the effect that leads to surface 

O • superconductivity, and it is not at first obvious why magnetic dots should support the relatively large supercurrents 

O ■ associated with multiple vortices. Therefore a possible oxide layer between the magnetic dot and the superconductor 

^ , may restore the boundary condition to the one of a superconductor/insulator interface. In the present paper we 

put a single magnetic dot on top of the superconductor and study the behavior of our sample in such non-uniform 

magnetic field of the dot which enhances the possibility of obtaining various combinations of superconducting states. 

i3 . To better understand the problem we start from a simple theoretical model for the inhomogeneous magnetic field 

" " ' profile that, we believe, captures allpaapects of the physics involved. Models used before vary from a representation of 

the uiagnetic dot by a perfect dipoldlB to a magnetic field profile calculated numerically for an infinitely thin magnetic 

diskEj. To obtain a better insight we start first with a simple step- like field model and subsequently investigate the 

more complicated real magnetic field profiles, which we obtained numerically. 

Theoretical studies have predicted that in mesoscopic disks surrounded by an insulating media three kinds of 
superconducting states can exist - giant vortex (a circular symmetric state with a fixed value of angular momentum), 
multi-vortex state (a collection of single vortices which can be .obtained as a linear combination of giant vortices 
with different angular momentum), and ring-shaped vortex state^3E2l with larger energy than giant and multi- vortex 
states. The ring-shaped two-dimensional vortex states have a cylindrically symmetric magnetic field profile and they 



are different from tfie ring-vortices which are e.g. found in three-dimensional supjerfluid hquid heHum. In the present 
paper we observe giant vortex states and first-order transitions between themu'cl, and for sufficiently large disks, 
multi-vortex structures, which are the analogue of the Abrikosov flux-line lattice in a bulk superconductor. The 
latter results not only from a mixture of giant vortex states but also from giant-ring vortex combinations as well. 
The latter one may even lead to an off-center location of a single vortex or multi-vortices. Moreover, with changing 
the strength of the field there is a transition between giant-giant and giant-ring multi- vortex states. Increasing the 
step height of the magnetic field profile we found re-entrant behavior, i.e., Iransition from giant to multi- vortex 
state and back to the giant vortex state before superconductivity is destroyedQ. We found that for sufficiently large 
magnetic disks vortex/ anti-vortex structures can be formed. In Qfder to investigate these different vortex structures 
we use the method proposed by Schweigert et alU and PalaciosQ'EI, with its semi-analytical extensions of Ref.EEl to 
determine the stability of the different multi-vortex configurations. In particular, the analysis of Rcf.E3 showed that 
in a superconducting disk the ring-shaped vortices are unstable in the presence of a homogeneous magnetic field.. 

Our analysis is within the framework of the phenomenological Ginzburg-Landau (GL) theory. Although this theory 
has only a firm mathematical derivation in a narrow range of magnetic field close to the superconducting-normal|State 
boundary, it has been found that it gives also very good results deep inside the superconducting phase diagranJj'El. 

The paper is organized as follows. In Sec. II we present our theoretical model. In Sec. Ill we discuss the giant 
vortex states and study the influence of a step magnetic field profile, with zero average which is centered at the disk or 
has a ring symmetry, on the superconducting state. These step profiles are limiting cases of the actual experimentally 
realizable profiles. The stability of multi-vortex states and transitions between them are investigated in Sec. IV. 
In the Appendix we give the analytical approach/solution to this subject. The excitement in this study is shown 
through different Hin — R phase diagrams in Sec. V, where i/i„ is the magnitude of the magnetic field profile and R 
represents the radius (of superconducting disk, magnetic dot, current loop, etc.). In Sec. VI we present the results 
of our approach applied to a superconducting disk in the presence of an experimentally realizable real magnetic dot 
and current loop field profile. Sec. VII is an extension of previous sections where the infiuence of an additional 
homogeneous background magnetic field is investigated. Our conclusions are given in Sec. VIII. 

II. THEORETICAL MODEL 

We consider a mesoscopic superconducting disk with radius R and thickness d surrounded by vacuum. The external 
magnetic field H = (0, 0, H) is directed normal to the disk plane. In this paper we investigate two different magnetic 
field profiles: 1) step-like magnetic field in the center of the disk, and 2) a ring step-like magnetic field profile with 
inner radius of the ring R^. The magnetic field strenghts of the profile (Fig. ||) are chosen such that the total 
magnetic flux equals zero. These models should correspond to the magnetic field of a perpendicular magnetized 
magnetic dot placed in the center of the disk on top of the superconductor, and the magnetic field due to a current 
loop placed on the superconductor, respectively. We assume that the superconductor will not alter the magnetic state 
of the ferromagnetic disk. We have to solve the system of two coupled non-linear GL equations which determine the 
distribution of both the superconducting order parameter, ^(T*), and the magnetic field (or vector potential A(T*)) 
inside and outside the superconductor 
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where the density of the superconducting current j is given by 

7 = £^ f^*^*-*^^*) -l^l^l^l". (3) 

im V / mc 

Here T* — {'p , z) is the three-dimensional position in space. Due to the circular symmetry of the disk we use 
cylindrical coordinates: p is the radial distance from the disk center, Lp is the azimuthal angle and the z-axis is taken 
perpendicular to the disk plane, where the disk lies between z = —d/2 and z = d/2. 
Eqs. (tt[]|) has to be supplemented by boundary conditions (BC) for ^(T*) and A{'1^): 



.i?i^_^:l]^ 



0, (4) 



where the subscript n denotes the component normal to the disk surface. The boundary condition for the vector 
potential has to be taken far away from the disk where the magnetic field becomes equal to the external field H 



A 



= 0. 



(5) 



Using dimensionless variables and the London gauge div A = we can rewrite the system of equations (|l|-y) and 
BC (0) in the following form 



(-iV -^] i/; = i/;- VlV'l^ 
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= 0. 
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Here all distances are measured in units of the coherence length ^ = h/ ^2m \a\, the order parameter in $o = vMT/^j 
the vector potential in c'h/2e^, k = A/^ is the GL parameter, and A = c-\/m/7r/4e^o is the London penetration 
depth. We measure the magnetic field in Hc2 = ch/2e£,^ = k^/2Hc, where He = y^Ana^/p is the thermodynamical 
critical field. 

The free energy of the superconducting state, measured in _fo = H'^V/8tt units, is determined by the expression 
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with the local magnetic field 



li {T') = "^ X ^(T"). 



We restrict ourselves to sufficiently thin disks such that d ^ ^, A. In this case, to a first approximation, the magnetic 
field due to the circulating superconducting currents may be neglected and the total magnetic field equals the external 
one Hq. Within this approximation we have to solve only the first GL equation (|g) with A = Aq, where Hq = rotAQ. 
One should notice that in this approach, for different cases, we change only the vector potential profile Aq. For our 
step profile 
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0, 0<p<Rd, 

Hin , Rd < P < Rl, 

— Hout, Rl < P < R2, 

0, R2< P<R- 



(10) 



the vector potential is given by 



^o(p) 



0, ^ < p < i?rf, 

^(p"^), Rd<P<Ri, 

IL^{Rl^Rl)-Ii^^p^f), R,<p<R2, 

0, R2<P<R- 



(11) 



where Hin describes the positive step and Hout = —HinR\/{Rj^ — R\) determines the value of the negative step (see 
Fig. 0, Fig. 1^ and the Appendix). 

First, we determine the z-dependence of "0 (^)- Representing the order parameter as a series over cosines '0 (T*) = 
00 
^ ipkilf) cos [kiTz/d) and using the same BC (H) at the disk sides {z — ±d/2) and using the first GL equation (Q), 

fe=0 

one can verify that the uniform part of the order parameter, i.e. the fc = term, gives the main contribution for 
(TT^/d)^ >> 1. Therefore, we may assume that the order parameter is uniform along the z direction of the disk 
and average the first GL equation over the disk thickness. After this averaging and for fixed value of the angular 



momenturn- it leads to ijj{j^) = f{p)e'Kjp{iLip), and the problem for / (p) is reduced to a one dimensional problem, 
like in Reffl: 

^ ^-P?f + {{--A] )/ = /(l-n, (12) 




and for the vector potential 



p dp dp 



i;^+S^) = (^^|^'«w«)'<^w/^'' ('^) 



where the function 6{x) = 1 (a; < 1), (a; > 1), and R, d are the dimensionless disk radius and thickness, respectively. 
The brackets (...) refer to averaging over the disk thickness. 

III. GIANT VORTEX STATES 

The giant vortex state has cylindrical symmetry and consequently the order parameter can be written as '0 {~f^) = 
f (p) exp (iLip). The stable states are obtained in the following way. From the linearized GL equation we find / (p) up 
to a multiplying constant. This function is then inserted into the free energy expression ffl) which after minimization 
determines: (i) the constant in /(p), and (ii) the energy value corresponding to-|the stable state. It can be shown 
that, for the case of giant vortex states, the present approach and the one of Reffl, which was based on a solution of 
the non-linear GL equation, result into the same functional F (Hin) dependence. 

The linearized GL equation for / (p) takes the form 

i/ = 0, L = ---^p^+(--Ao) -I. (14) 

pdp dp \p J 

The superconducting state starts to develop when the minimal eigenvalue of the operator L becomes negative. For the 
zero angular momentum state, the normal state transforms to the superconducting one with decreasing magnetic field 
below the nucleation field Hnuc- For nonzero angular momentum, the superconducting state appears when we cross 
either the lower Hnuc,i or the upper Hnuc,u critical magnetic field which depends on the disk radius. The eigenvalues 
and eigenfunctions of the L operator are found from 

Lfn,L (p) = K,Lfn,L (p) , (15) 

where fL,n{p) satisfies p{df/dp)\Q = at the disk center. The index n = 1,2,... enumerates the different states 
for the same L-value. 

In general, the eigenfunctions of Eq. ( [l5[ ) can be obtained analytically in the case of our step magnetic field profile. 
We present the complete calculation in the Appendix for both considered cases. Alternatively, we solve Eq. ( [l5| ) 
numerically through the finite difference technique. We put the order parameter on a space grid and find numerically 
the eigenfunctions and eigenvalues of the operator L using the Householder technique. 

We start our analysis with the magnetic field profile shown in Fig. |l|. First, we consider a magnetic dot on top of 
the center of the disk, i.e. with R^j = 0, and R2 = R — 6.0^ and study the influence of i?i, i.e. the width of the 
positive field region, on the superconducting state. 

The magnetic field dependence of A for different angular momenta L are shown in Figs. 0(a-c) for the lowest radial 
state, i.e. n — 0, and in Figs. 0(d-f) for the first radial state, i.e. n — 1 for three different values of i?i. The top 
axis shows the fiux corresponding to the positive magnetic field region 0i„ = HmTrRi, which in dimensionless units 
becomes 4>inl4>a — {Hin/2Hc2){Ri / £,)'^ where 00 — ch/2e is the quantum of flux. All numerical calculations were done 
for a disk thickness d/^ = 0.1 which is within the thin disk approximation. The negative L values (dashed curves 
in Fig. B) correspond to "anti- vortices" in conventional superconductors. The thin horizontal line gives the A = 
level. From Fig. g one notices that with increasing i?i: 1) the eigenvalues A of the states with the same L become 
more negative; 2) the magnetic field range over which solutions of Eq. (nj) can be found decreases; 3) the number of 
possible solutions decreases also; and 4) the A(iJi„) dependences become more parabolic. The latter can be explained 
by the fact that increasing £i_rorresponds to a more homogeneous magnetic field profile inside the disk, i.e. we reach 
the case considered in Refs.l3o. 

For small i?i the curves L and L + 1 anti-cross for sufficiently large L values and consequently the low vorticity 
states have lower energy even with increasing strength of the magnetic field profile. E.g. for Ri/S, = 1.5 this occurs 



when L > 4 = L* and for Ri/$, = 3.0 when L > 5 — L* . Notice that the slope of A for the L > L* curves is 
substantiaUy smaller than for the L < L* curves in the high field region. Notice also that the n — 1 states have a 
higher A value than the n = states and their energy is also larger than those of the anti- vortex states for the same 
vorticity. 

In Figs. 0(a-c) the radial dependence of the superconducting density j-^l is shown for \L\ < 3 at Hin — 0.75i?c2 
for the corresponding profiles of Fig. g. Notice that, with increasing Ri, the Cooper-pair density near the edge of 
the sample becomes more non-homogeneous. The \L\ ^ states have a vortex sitting in the center of the disk, i.e. 
iIj{p = 0) = 0, which becomes larger with increasing \L\. For the situation of Fig. ^(c) there is a narrow, very negative 
magnetic field region of Hout/Hc2 = 0.9643 in the region 4.5 < p/£, < 6.0 which leads to a considerable suppression of 
the Cooper pair density. For example, the L = vortex state, i.e. the Meissner state, has a strongly reduced Cooper 
pair density for p/$^ > 1. For the excited state {n,L) = (1,0) (dashed cupips in Fig. ^) the order parameter vanishes 
inside the disk and a ring-shaped node in the wave function ip is formedE2lEil which leads to a ring-like vortex. 

The eigenvalues A determine the free energy F of the giant vortex state. For the giant vortex state we consider 
only states which lie below the F = level. In this approximation the order parameter is 
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and the minimal energy value is 



where 



^{lf)=(-A^-^) /„x(p)exp(zL^), (16) 



F = -A' — ^, (17) 
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h= j pdp ft^L (P) > ^2 = / pdp fl^L (P) ■ 



The dependence of the free energy on the magnetic field strength of the inner core of the magnetic field profile, Hin, 
are shown in Figs. §(a-c) for the angular momenta \L\ < 11 for different values of i?i. The free energy is expressed in 
units of Fq = H'^V/8n. The highest value of vorticity in this disk is L = 15 (see Figs. p|(a-c)). From a comparison 
of the magnetic field dependence of F (Hin) for Ri/S, — 1.5 (Fig. H(a)) and the one with Ri/S, — 4.5 (Fig. H(c)) we 
clearly observe the reduction of superconductivity and the reduction of the maximal possible vorticity with increasing 
i?i . The envelope of the lowest parts of the curves in Fig. H represents the field dependence of the ground state 
energy. Notice that the increase of the width of the positive magnetic field region leads to an increase of the energy 
of the ground state. With increasing applied field the L — > L + 1 vortex transitions take place at the field where 
the corresponding curves cross (for example, the — > 1 transition occurs at Hin — 0.8705i/c2 for -Ri/C — 1-5). The 
crossing points are shifted towards lower field values when increasing Ri . The L — > L + 1 transitions are of first order 
and lead to jumps in the magnetization of the sample. Notice that the positive flux captured in the superconducting 
disk for different L-states is not quantized. This is made very clear in the (pin — Ri diagram presented in Fig. H(d) 
which shows the ground state vortex configurations. Notice that the flux in the positive magnetic field region has 
to increase with more than one flux quantum (Jjq before the vorticity of the superconducting state increases with cpe 
unit. A similar phenomena was observed earlier for mesoscopic disks and rings in a homogeneous magnetic fieldE3. 
Notice also that for higher L values quantization is slowly restored, and Ac^ decreases with enlarging i?i, i.e. the 
radius of the positive field region. 

Following our assumption that field,inside the superconductor equals the external one, we obtain the expression for 
the superconducting current densit}£j: 

J„x-%^(--Ao)/,U(p) (18) 

An,L P 

where A„_i is determined by Eq. (p^), Aq represents the vector potential of the applied field and 

R R 
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An,L ^ — pdp fn^L (p) , Bn,L ^ — pdp /„^i (p) . 


The magnetic field due to the supercurrents, neglected in our first order approximation, is calculated from 



rotiC = -^7- (19) 

Since the sup ercu rrent has only an aziniuthal component, and, is situated only in the superconductor plane, the p 
component of Hsc can be neglected. Consequently, we obtain 

HUp) ^ --^ J jMdp, (20) 

and the magnetization of the superconductor is then defined as the magnetic field expelled from the superconductor. 

M = J El2lfl_I^dV = \ J HUp)pdp, (21) 



where Hq denotes the applied magnetic field, and d is the disk thickness. 

The corresponding M {Hin) curves are given as insets in Figs. E|(a-c). The phase transition from the superconducting 
to the normal state is of second order (all curves F {Hin) reach the F — line with zero derivative). The curves 
F {Hin) in Fig. || which are situated above the ground state energy correspond to metastable giant vortex states. 
With increasing applied field the transition from the Meissner state {L = 0) to the normal state goes through a set of 
consecutive first order transitions between the L and L + 1 giant vortices which is finished by a second order transition 
to the normal state. 

Our next step was to fix the magnetic field profile (we took Rd/S, = 0.0, Ri/£, = 4.5, R2/£, = 6.0) and enlarge 
the disk radius R. In Figs, ^(a-b) the radial dependence of the superconducting density \iIj\ is shown for \L\ < 3 at 
Hin = 0.75Hc2 and for a superconducting disk of radius R/(, = 9.0, and 12.0, respectively. One can see that, even for 
L = 0, superconductivity in the center of the disk is destroyed, which is opposite to the homogeneous magnetic field 
case where |^| is maximal at p = OB. The shape of the L = curve changes drastically and this occurs already for 
small enlargement of the disk. This transition is shown in Fig. ^(a). By increasing R/^ from 6.1 to 6.5 the Cooper 
pair density in the center of the superconducting disk decreases from 0.73 to 0.03. Furthermore, in Fig. 0(b) we see 
also qualitative changes in the Cooper pair density of the excited state {n,L) = (1, 0) (dashed curve). Now the order 
parameter vanishes twice inside the disk and a double ring-like vortex is formed. This transition is shown in Fig. Hb) 
with increasing disk size from R/^ = 10.0 to 11.0. 

The dependences of the free energy on the magnetic field Hin for different sizes of the superconducting disk are 
shown in Figs. CTa-d). Notice that there are large differences in comparison with the previous cases. We still have 
the L — > L + I transitions, which are of first order and which lead to jumps in the magnetization of the sample (see 
Figs. @(a-c)), but with increasing applied field the transition from the Meissner state {L ~ 0) to the normal state 
goes through a set of consecutive first order transitions and is finished by a first order transition back to the L = 
state, sometimes with an intermediate L = 1 state. For R/^ = 9.0 we have 0^-1^2^-3^1^-0 transitions, for 
R/£, — 12.0 we have 0— >1— >2— >l-^0, and for R/^ = 18.0 only ^- 1 ^- 0. Moreover, for a sufficiently large disk 
(see Fig. M(d) which corresponds practically to the R/S, — > oo situation) we have no first order transitions but only 
the Meissner state as the ground state. The reason is that for i?/^ -^ oo the non-zero magnetic field region is limited 
(relatively) to a small area and we can always define a circle with radius p sufficiently large where the superconducting 
current is zero and, consequently one must have L = 0. But, this Meissner state is qualitatively different from the 
"usual" one. The radial distribution of the Cooper pair density is extremely inhomogeneous and superconductivity is 
strongly suppressed in the interior of the disk (see Figs. H(b) and @(a)). 

IV. MULTI- VORTEX STATES 

It is well known that for sufficiently large disks the giant vortex state can break up into multi-vorticesLTEl. In our 
analysis, we considered different field profiles and disk geometries to investigate the properties of this transition. For 
smaller disks, the confinement effect dominates and we found that only the giant vortex states are stable and possible 
multi-vortex states, if they exist, have alwjays larger energi^. In order to investigate such structures in our case we use 
the method proposed by Schweigert et alU andrEalacios BB and estejad|ed it to determine the stability of the different 
multi-vortex configurations as proposed in Ref.EJ. Following Refs.ElBcj the order parameter of the multi-vortex state 
is written as a linear combination of eigenfunctions of the linearized GL equation ( [14|) 

L oo 
L,=On=0 



where L is, in the homogeneous magnetic field case, the value of the effective total angular momentum which is equal 
to the number of vortices in the disk, and n enumerates the different radial states for the same Lj . Later on, we will 
see that for our inhomogeneous magnetic field case the assignment of the total vorticity can be more tricky. 

Substituting ( |22| ) in the free energy expression (|9|) we obtain i^ as a function of the complex parameters {Cn,Lj}- 
Minimization of F with respect to these parameters allows us, to find the equilibrium vortex configurations, and 
to determine their stability. We use the procedure described in Rcfs.EJt3 with full consideration of combinations of 
vortex configurations with different radial states. 

We begin our analysis with states built up by only two components in Eq. (|2^). This brings quantitative bounds 
in our calculation but, nevertheless, will give the correct qualitative behavior and facilitates the physical insight into 
the problem. The free energy of a two component state built out of (ni, Li) and (712, ^2) is 

F — ^m.Li^niXl + ^n2,L2^n2,L2 + 4C„^ j^^ C^^^Lj^m ,n2,il ,-^2 + '^^ni,LiC„^i^^Bn^^L^ + 2A„2_i2 C'„2,L2 -^"2X2 ; (23) 

where 
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Ah,L, = — / Pdp fn^x, (P) ) 

R 



An,,n2.Li,L2 = — / pdp fl^^L^ (p) /n2,i2 (P) ^ 
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Bn,,L, = -^ I pdp fn^^L, (p) • 
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One should notice that we leave the possibility of combination of states with different radial states, i.e. ni and n2. 
Although, in general. Cm, Li is a complex number, for our two component state Cm, Li is a real number. Minimization 
of Eq. (p3) with respect to C^.l^ and Cn2,L2 gives for the multi-vortex states: 

/ \ ^^^ 

^(0) if ^J^ni,LiAn2,L2Bnx,Lx + ^■'in2,L2^ni ,n2,ii,i2 ^"2 ,^2 \ (c^.\ 

y Ani,LiAn2,L2 *^ni,n2,Li,L2 / 

/ \ 1/2 

^(0) I / ^■''-n2,L2^-£'l, "1-^12,-^2 "I" ^^^Li,niAni,n2.Li,L2^Li,ni \ 

^"2,^2 ~ 4 7^ 'r _ Zl 42 I ' 

Y ^ni,Li^n2,L2 ^^ni.n2,Li,L2 ) 

and inserting these expressions into Eq. (E3) leads to the energy of the multi-vortex state 

jp __ ~^rii,Li^"2,i2-°ni.Li ^ ^^2,^2 "I'-^i "2,^2 '^ "I'-^i "2,i2^ni:n2,il,i2-°ni,Li-D„2,L2 /„p_N 

^ni,«2,Li,L2 - -1 —7 -4 42 '-^°'' 

^ni,Li-^n2,L2 ^^ni.n2,Li,L2 

The corresponding conditions for the stability of the vortex state are 

O r o^ni,Li v ^^n\,L\-^n2.L2^n-i,L\ ' ■^^^n2,L2^ni,n2,Li,L2-^n2;L2 



'^'^'ni,Li Am,LiAn2,L2 - ^Aii,n2,Li,L2 

Or oJ\-yi2,L2 V -^^n2,L2 ^ni ,Li -^712 ,L2 ' ^-^^ni,Li^ni,n2,Li,L2-^ni,Li ^ 



> 0, (26) 

>0, 



^^n2,L2 An^,LiAn^,L2 ^ ^Aii,n2,Li,L2 

d^F d^F ( d'^F y 

-""2,^2 / 
r A.^r r n. r \ 

{~-^n2,L2Ani,LiBn^,L2 + 2A„j ^^i ^m ,n2,-tl ,^2 -^"l ,-tl ) > 



^'^ni,Li ^'^n2,L2 \dCm,LidCn2,L 

64 ( — A„j^Li^n2,i2^ni,il + '^^n2,L2A-ni,n2,Li,L2Bn2,L 



A T A r — 442 
^ni,Li^n2,L2 ^^m ,n2,-Li ,-L2 

In our analysis we investigated the influence of the width of the positive field region and radius of the disk on the 
phase diagram, and especially on the stability of the multi- vortex states. For small values of Rxj^ these multi- vortex 



states are always metastable. However, with enlarging Ri/S, these states can obtain lower energy. The energies of 
the equilibrium vortex states are plotted in Fig. ||(a) for Ri/S, = 5.25 and i?2/C = -^/C = 6.0 and in Fig. g(b) for 
i?i/^ = 4.5 and i?2/C = 6.0 with a larger disk radius R/£^ = 9.0. The giant vortex states are given by solid curves, anti- 
vortex states by dashed and the multi-vortex states (dotted curves) by (Li : ^2), i-e. the angular momentum values 
they are composed of. As expected, when enlarging the disk, multi-vortex states become more stable and, moreover 
in Fig. y(b), we observe the existence of multi- vortices as a combination of a giant and a ring- like vortex. These 
states have surprisingly low energy and become the ground state for specific ranges of magnetic field. Furthermore, 
when increasing the field, we observe a phase transition between this giant-ring and giant-giant multi-vortex states, 
when another type of multi-vortices becomes the ground state. The giant-ring multi- vortex states with lowest energy 
are given in Fig. ^(b) by the quantum numbers ((ni,Li) : (712,^2)). It should be noted that there are many other 
metastable combinations possible, which are not shown in the figures. 

As shown in Fig. 0, with increasing i?i/^ and the size of the superconducting disk the multi- vortices become more 
stable and they can even become the minimum of the F {Cni,Lx, ^^712,^2) function. We focus our attention to Fig. pKa). 
The solid and dashed curves represent the giant vortex and anti- vortex states, respectively. The dotted curves are the 
energies of the multi- vortices. For example, lets follow the L — curve which splits into the (0 : 4) multi- vortex state, 
at Hin/Hc2 — 0.28 and with further increase of the field becomes another giant vortex state, but now with L = 4, at 
Hin/Hc2 — 1-0. The concomitant change of the Cooper pair density is illustrated in Fig. |l^. Vortices enter the disk 
from the boundary and move to the middle with increasing field, and join into a giant vortex again. 

For a fixed magnetic field profile and with increasing size of the superconducting disk, we obtain a variety of different 
superconducting states. Let us discuss the giant-ring configurations first. We observe transitions between the giant 
and multi- vortex states with the same vorticity. The contour plots of the |V'| distribution for different multi- vortex 
states as ground states (((0, 1) : (1,0)), ((0,2) : (1,0)), ((0,3) : (1,0))) are shown in Figs. |ll|(a-c) for a magnetic field 
profile with Ri/^ = 4.5, i?2/C = 6.0 and i?/^ = 9.0, for different values of Hin (dark regions correspond to high density 
and white regions to low density). The most interesting result we obtain for the ((0, 1) : (1,0)) combination. This 
configuration becomes the ground state configuration for Hin/Hc2 — 0.75. Surprisingly, the contour plot of the \ip\ 
distribution is not circular symmetric - there is one vortex positioned off-center. In the lower part of the same figures 
we give the corresponding contour plots of the phase of the superconducting wave functions. The contour plot of the 
phase also nicely illustrates the off-center location of the vortex. For the two other cases more vortices are present 
which are located on a ring centered around the center of the superconducting disk. The giant-ring multi- vortex state 
can be seen as a transition between (1, 0) and (0, 5) state (see Fig. H(b)). As function of the magnetic field we start 
from a ring vortex, for Hin/Hc2 = 0.03, this state splits into a giant-ring vortex state, and finally, for Hm/ Hc2 = 0.98 
we obtain a giant vortex state. This remarkable phenomenon is illustrated in Fig. O. 

However, giant-giant vortex combinations show completely different behavior. As shown in Fig. ^(b) with increasing 
magnetic field i?m, combinations (0 : 7), (0 : 8), (0 : 9) etc. become the ground state. In the contour plot of the |'i/'| 
distribution we observe single vortices which are arranged on a ring with a low density area situated in the center 
of the disk (see Figs. |3|(a-c)). This central area is associated with a giant vortex, and encircling it leads to a phase 
change showing the vorticity 7. Fig. 0(d) shows the contour plot of the corresponding phase with seven anti-vortices 
arranged in a circle around the giant vortex and the total vorticity is L = 0. Figs, ^(e) and Fig. |l3|(f) correspond, 
respectively, to a giant vortex state with vorticity 8 and 9 where, respectively, 8 and 9 anti- vortices are located between 
the giant vortex and the boundary. Notice that the Cooper pair density for the L = state for this disk geometry 
(i?i/^ = 4.5, -R2/C = 6.0, R/^ = 9.0) shows a similar behavior with a low density area in the center of the disk 
(see Fig. ||(b)). It is obvious that enlarging the superconducting disk enhances the influence of the negative part of 
the step-like magnetic field. In the case when multi anti-vortices are involved, the total vorticity of the giant-giant 
multi- vortex state is equal to the lowest vorticity of the two giant vortex states which compose the vortex state. 

One more result should be noted. Following the ground state free energy diagram (Fig. |^(b)) we notice that the 
total vorticity doesn't change uniformly (0^1 -p^2 — >3^1^7— >8— >9) and that the jumps in vorticity |AL| 
are not always equal to 1. This differs from Ref. cil where it was claimed that the lowest barriers are those between 
the L and L ± 1 states. We found that the L — > L ± 1 transitions take place in disks for a small maximum value 
of vorticity or at magnetic fields close to the "superconducting-normal state" transition point. Between these two 
limitirig regimes L — > L ± N transitions are possible with |AL| = A^ > 1. Our results are also in agreement with 
Rcf.cj where it was found numerically that several vortices can enter (or exit) at once for disks with sufficiently large 
radius. In Sec. Ill, where we investigated the inffuence of Ri/^ (with i?2/C = R/S, — 6.0) we only found L — *■ L + 1 
transitions. 

In order to show the complexity of the system under study, we investigated one more field profile, the one shifted 
from the center of the disk i.e. a ring magnetic field. We keep all parameters from the previous case, and shift the 
field by Rd — 2.0^ (i?i — Rd ^ 4.5^, R2 — Rd = 6.0^, R/S, — 9.0) towards the disk edge (a top view of this profile is 
given schematically in the inset of Fig. O). As shown in Fig. O, with Rd ^ the multi- vortices become more stable 



and for Hin/Hc2 > 0.31 are the ground state. The sohd curves represent the giant vortex states, dashed curves denote 
the energy of the anti-vortex states while the dotted curves correspond to the energy of the multi-vortices. Since 
giant vortex states show similar, re-entrant behavior as for the Rd = case, the multi-vortex configurations behave 
analogously. Giant-giant multi-vortex states dominate the free energy diagram while the giant-ring combinations are 
present only as metastable states and are not shown in Fig. n3. However, one should notice the presence of giant-multi 
anti-vortex states, again strongly correlated with re-entrant behavior. The complete equilibrium phase diagram for 
Rd G [0.3, 3] is given in the next section. 

We checked that for the disk parameters which we used, an increase of the number of components in Eq. ( |22| ) does 
not lead to different vortex configurations in the ground state. In order to investigate the region of stability of multi- 
vortex states in the above analysis we took the order parameter as a linear combination of three components in Eq. ( |22| ) 
and minimized the free energy with respect to the three variational parameters C^. . The giant and multi- vortex states 
considered before correspond to the extrcmum points of the F{Cli j ^L2 j C'lj) function. This analysis has shown that 



accounting of third component in Eq. (22) gives all states obtained in two component consideration. Additionally, 
it results in 1) possible reducing of the region of existence of metastable multi- vortex states at the low magnetic 
field limit, 2) appearance of additional unstable states corresponding to the sadle points of the F(Clj, Cl^, Cig) 
function, and, 3) no new vortex configurations in the ground state. We did a similar investigation using a numerical 
approach of Schweigert and Peeter£3. This was done for five component vortex configurations, and same results were 
obtained. Using two different approaches, for three and five components, no new ground state configurations are 
found. Moreover, energies of same states found in both analysis for different number of components differ less than 
0.2%. Because all multi-vortex configurations considered in this paper are in the ground state, or near it, we conclude 
that two components in Eq. (E3) are enough to describe the vortex structure in our system. 

V. HiN - R PHASE DIAGRAMS 

First, we investigate the influence of the width of the positive magnetic field region on the different vortex configu- 
rations. Having the free energies of the different giant vortex configurations for several values of Ri/^, we construct an 
equilibrium vortex phase diagram. Fig. |lj shows this phase diagram for a superconducting disk with radius R = 6.0^ 
and thickness d — 0.1^ where we took i?2 — R- The dashed curves indicate where the ground state of the free energy 
changes from one L state to another and the solid curve gives the normal/superconducting transition. Notice that 
the superconducting/normal transition moves towards lower fields with increasing radius of the positive field region. 
Also notice that we don't have any negative L (anti-vortex) state as ground state even when the negative field area 
is much larger than the positive one. 

As shown in previous section, an increase of Ri/£, is able to bring the energy of the multi- vortex states below those 
of the giant vortex states and they can become the ground state (see Fig. ^(a)). In our phase diagram, the area 
bounded by the thick curves denotes the region of existence of the giant-giant multi-vortex states as ground state. 
One can see that, with enlargement of the positive field region, ground state multi-vortices appear for i?i/^ = 4.42 
and Hin/Hc2 = 0.91. Further increase of Ri/£, broadens this multi-vortex area to almost the whole superconducting 
region. These multi-vortex states consist of vortices in a ring structure, and the total vorticity equals the highest 
vorticity of the giant vortex states involved (at least in the case when no anti- vortices are involved). 

In order to show that the stabilization of the multi-vortex states due to an inhomogeneous magnetic field is not 
peculiar to the R = 6.0^ disks, we repeated the previous calculations for a larger superconducting disk. Furthermore, 
we investigated in detail the phenomena shown in previous section (Figs. ^, |l2|, |l3|) . 

The effect of the size of the superconducting disk on the phase diagram is illustrated in Fig. Qq. The parameters 
considered are Rd/^ = 0.0, i?i/C — 4.5, and i?2/C = 6-0. The solid lines indicate where the ground state of the free 
energy changes from one state to another (either giant or multi- vortex state) and dashed lines correspond to transitions 
between different giant vortex states as metastable states. One can clearly see the re-entrant behavior. For example, 
for i?/^ = 9.0 we observe the change of the total vorticity asL = 0— >1-^2^3^1^0. The L = and L = 1 
states as ground state cover the largest part of the phase diagram. With increasing disk size, all other giant vortex 
states are strongly suppressed in favor of the various multi-vortex states. Different giant vortex states are present as 
the ground state for small disk sizes but with increase of disk size, islands with different multi-vortex configurations 
{giant-ring and giant-giant multi-vortex states) dominate the ground state diagram. Precisely, for R ~ 6.38^, giant- 
ring multi-vortex state appear as the ground state, and for R — 6.77^ we obtain the giant-giant multi-vortex states 
as the ground state. In the previous diagram in Fig. |lj we have shown the existence of the giant-giant multi-vortex 
states as the ground state for R = 6.0^. With slight increase of the disk size, these states become metastable, and 
with further enlargement of the superconducting disk, they become the ground state again. However, these states 
are different from the previous ones and in the contour plots of the Cooper pair density of these "new" giant-giant 



multi-vortex states (see Fig. |3|) we found a giant vortex state in the middle of the superconducting disk surrounded 
by anti- vortices, and where the total vorticity is now equal to the lowest vorticity of the giant vortex states that the 
multi-vortex consists of. For R > 12.33^, the multi-vortex states become metastable again and the L = and L — 1 
states become the only ground states. Moreover, for R > 19.34^ we have the Meissner state for all values of the 
applied magnetic field. 

To present the complexity and excitement in this study, we give one more diagram. As shown before, moving the 
step-like field profile along the radius of the superconducting disk, i.e. a ring magnetic field profile as in the case of 
a current loop, stabilizes the giant-giant multi-vortex states as the ground state. In Fig. [l^, we present the phase 
diagram as function of Rd- The parameters of the magnetic field profile are i?i — Rj. = 4.5^, i?2 -^ Rd — 6.0^, and 
R/^ = 9.0. Thick solid curves indicate the transitions between different vortex states and dashed lines denote the 
transitions between metastable giant vortex states. The re-entrant behavior, as in previous case, is clearly visible. 
The most important result is that the multi- vortex states dominate this diagram. Moreover, we have both types of 
multi- vortices, i.e. giant-giant and giant-ring, as stable and ground state, and with shifting the field profile towards 
the disk periphery, giant-giant multi- vortex configurations cover most of the superconducting region. . The giant-giant 
multi-vortex states in this case appear both as states with no anti-vortices present, and as giant-multi anti-vortex 
states. Following the transition lines, a correspondence between different giant-multi anti-vortex states (Li:L2) can 
be seen: L2 remains the same, while Li increases from to 2. As one can see, the latter is strongly correlated 
with re-entrant behavior. In the rest of the diagram, the multi-vortex states with the "classical" geometry (both 
giant-giant and giant-ring) dominate. However, a difference between those two states exists. As shown in Fig. y_2|, in 
the case of a giant-ring state, vortices are preferentially distributed within a ring-shaped lower density area. Further, 
considering the high density areas at the disk periphery, a shift in phase of A9 = tt/L is observed, in comparison 
with corresponding giant-giant state (see, for example. Fig. p^ ), where L is the total vorticity. Although these states 
have a different origin, sometimes they can exhibit a similar distribution of vortices, but the phase is always able to 
distinguish between them. 

VI. COMPARISON WITH REAL MAGNETIC FIELD PROFILE 

As emphasized before, our step-like field model is a simplification of the magnetic field profile of a ferromagnetic 
dot (see solid curve in Fig. |l]) as shown schematically in Fig. |8[ In this case, the radius of the dot corresponds to i?i, 
and i?2 is always equal to R. The magnetic field and vector potential were calculated through numerical integration 
of 

-Ri 



H,{p,z) = m, / dr' , ^'^'^^^ ^^ E 

' ^(r' - p)2 -I- (z' _ z)2((r' + p)2 + (z' - ^)2) 




and 



Ar'p 



(r' -|-p)2 + (z'- z)2 



(27) 



p 
A^{p,z) = - / rH^{r,z)dr, 


where ruz denotes the magnetic moment of a dot directed along the z-axis, i?i is the radius of the dot, Zq gives the 
distance from the plane of interest, and dd = z[ — z'q is the thickness of the dot. E{x) is the complete elliptic integral 
of the second kind. The calculations were done for z ^ 0, which allowed us to neglect the radial component of the 
magnetic field. 

In order to show the quality of our previous model we repeated the analysis from previous sections. First, we enlarge 
Ri, i.e. the radius of the magnetic dot, investigating the influence of the magnetic field on the vortex structure and 
especially on the stability of the multi- vortex states, and secondly we studied the vortex configurations resulting from 
an increase of the superconducting disk size. 

For small values of Ri/^ the multi-vortex states are always metastable states. However, with enlarging this param- 
eter these states lower in energy. The energies of the equilibrium vortex states as function of magnetic moment of the 
dot measured in units of toq = -ffc2, are plotted in Fig. |9|(a) for radius of the dot i?i/^ = 4.0 and R/S, = 6.0 and in 
Fig. |l9|(b) for Ri/^ — 4.0 with a larger disk radius R/^ = 9.0. Dashed curves correspond to anti- vortex states and dot- 
ted curves represent the energy of the multi- vortex states. The giant-giant multi- vortex states are given by (Li : L2), 
i.e. the angular momentum values they are composed of, and the giant-ring multi-vortex states with lowest energy 
are given in Fig. ^^(b) by {{ni,Li) : (ri2, £2))- It should be noted that there are many other metastable combinations 
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possible, which are not shown in the figures. Notice from Fig. |9|(a) that the results are qualitative similar to the 
results obtained with our step magnetic field model (see Fig. g). The difference is caused by the fact that part of the 
negative magnetic flux does not penetrate the superconducting disk. Because the positive field region in the center 
dictates the behavior of the phase diagram even when the total flux is zero (see, for example. Fig. |^(a)), we obtain 
very similar results with the real magnetic field profile, but with an increased number of possible superconducting 
states. However, increasing the size of the superconductor brings some qualitative changes. First, it is clear from 
Fig. &y(b) that there is no re-entrant behavior, and, second, no giant-giant multi-vortex configurations are ground 
state. But, this was expected since in our model field profile the whole magnetic flux is trapped in the center of the 
disk, with total flux equal zero. For the real profile, the total flux would be zero if our disk is infinitely extended, and, 
more importantly, the flux is now spread over the whole disk area. Naturally, smaller magnetic dots in combination 
with large superconductors would make the correspondence better, since most of the flux would be captured inside 
the disk. 

A similar discussion holds for the results for the real magnetic field profile of a current loop placed on top of a 
superconductor. This field profile (solid curve) is compared with our model in Fig. |2y (thick dashed curve). The 
magnetic field and vector potential were calculated numerically from 



H.{p,z) 



Ik 



r2 J2 



(i?l - p)2 + Z- 



:E{k^)+K{k'') 



and 



A^{p,z) = 




K{k^)-Eik'^) 



(28) 



(29) 



with 



k = 2. 



[Ri 



RiP 



■ /9)2 + Z2 ' 



(30) 



where K{x) is the complete elliptic integral of the first kind and i?i denotes the radius of the loop with current /. 

The free energy is shown in Fig. ^ for a loop with radius Ri — 7.5^, with superconductor disk size R — 9.0^, as a 
function of current / measured in units of /q = 7r^i?c2/Mo- Typical values of /q are 3.29mA, for aluminum, to 0.823A 
for high-temperature superconductors. One should compare these results with those presented in the Hin — Rd phase 
diagram (Fig. UTi) and we find that our previous model contained all the essential physics of the system. The free 
energy diagrams of both Fig. n^ and Fig. E^ show the re-entrant behavior in total vorticity and the existence of the 
giant-giant multi-vortex configurations over a large region of the phase diagram. 

VII. SUPERCONDUCTING DISK WITH A MAGNETIC DOT OR CURRENT LOOP ON TOP OF IT IN 
THE PRESENCE OF A BACKGROUND HOMOGENEOUS EXTERNAL MAGNETIC FIELD 

In addition, we investigated the vortex structure of a superconducting disk in the presence of an inhomogeneous 
magnetic field profile resulting from a ferromagnetic dot and a homogeneous external background magnetic field. In 
Figs. E2(a-c) we present the free energy as function of the external homogeneous field, for different magnetic dot 
thicknesses. The parameters were - radius of the dot i?i/C = 4.0, radius of the superconducting disk R/^ = 6.0 
with thickness d = 0.1^ and fixed ferromagnetic dot profile which is shown as inset of Figs, pl(a-c). As expected, 
when a negative external field overwhelms the average of the positive magnetic field of the dot in the disk center, the 
anti-vortex states become energetically more favorable, and, opposite, when the positive external field becomes larger 
than the average negative value of the ferromagnetic dot field, we see that the ground state goes through successive 
giant vortex states towards the normal state. For example, for a thickness of the magnetic dot dd — O.lf we obtain 
an almost symmetrical figure with respect to Hext — — 0.14iJc2- However, with increasing magnetic dot thickness, the 
magnetic field of the ferromagnetic dot (given as insets in Figs. p2|(a-c)) becomes more pronounced and we obtain two 
sets of curves corresponding with the vortex and anti-vortex states, each of which having two minima. These local 
minima occur at an external field value which is approximately equal to the average field of the positive or negative 
region of the magnetic dot profile. The region between these two minima is characterized by a strong interplay of 
states. 

For a current loop on top of the superconducting disk in a homogeneous external field we obtained qualitatively 
similar results. Also two sets of curves, with two local minima are visible in the free energy, for sufficiently large 
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currents in the loop (see Figs. ^(b-c)). We show numerical results for a system consisting of a superconducting 
disk with i?/^ = 9.0 and thickness d = 0.1^, and current loop with radius i?i = 7.5^ (see Fig. gO]). Just like in 
previous case, minima in the free energy are related to average fields in positive and negative region of the current 
loop magnetic field profile. However, several differences exist with previous case. First, with increase of the thickness 
of the magnetic dot in previous situation, the magnetic field profile becomes more similar to our step-like field model. 
For the current loop case this is not so, and with increasing current in the loop the magnetic field becomes strongly 
inhomogeneous, rapidly increasing in the vicinity of the loop. As a consequence the minimum in the free energy for 
negative external field is lower in energy, and, moreover, the other minimum slowly diminishes with further increase 
of the current in the loop. In this case even for negative applied external field, the positive peak in the magnetic 
field profile of the current loop is responsible for the stabilization of the positive L-states, up to large Hext < (see 
Fig. p3|(d)). It is obvious that the interplay between the external and the current loop magnetic field exhibits more 
interesting physics due to the strong field inhomogeneity. 

VIII. CONCLUSION 

We studied the superconducting state of a thin superconducting disk in the presence of an inhomogeneous magnetic 
field. A "model" step-like magnetic field profile was considered which is an approximation for the magnetic field profile 
resulting from a magnetic disk or from a current loop. The superconducting disk is assumed sufficiently thin that the 
magnetic field produced by the superconducting currents can be neglected. The effect of the width of the positive field 
region, the position of the field and the size of the superconducting disk on the vortex configuration was investigated. 
Numerous phase transitions were found, between states with different angular momentum number and between giant 
and multi-vortex states, plus transitions of the ground state between different multi-vortex configurations (giant-ring 
and giant-giant multi- vortex states). The model step- like magnetic field (total flux through the system equals zero) 
is found to stabilize the multi- vortex states both as the ground state (i.e., with minimal energy) and as metastable 
states. Increase of the width of the positive field region enhances the stability of the giant-giant multi- vortices, while 
enlarging the superconducting disk decreases the energy of the giant-ring multi-vortex state. In this case, the giant- 
giant multi- vortex states can also be the ground state but it represents a new configuration - a giant vortex surrounded 
by anti-vortices, and the total vorticity equals the lowest vorticity of the giant vortex states which are involved. We 
found that with an increase of the disk size, re-entrant behavior of the total vorticity is possible. This behavior can 
also be seen by shifting the field profile towards the disk boundary, i.e. for the current loop magnetic field, when 
giant-giant multi-vortex configurations cover most of the superconducting phase diagram region. Comparison with 
results obtained for a real magnetic field profile shows very good correspondence with our step-like models for the cases 
when most of the magnetic fiux is captured inside the superconducting disk. For larger radius of magnetic dot, the 
positive magnetic field region of the inhomogeneous magnetic field profile dominates and with increasing the magnetic 
dot radius, the total flux in the superconductor also increases, which is different from our model where the total flux 
always equals zero. Thus, the positive part of the magnetic field profile determines mainly the superconducting states, 
which leads to slightly different physics - no re-entrant behavior is present, more different superconducting states are 
possible, but without the giant-giant multi-vortex states with anti-vortices, since they are strongly correlated with a 
re-entrance of the vorticity. 

Addition of the background homogeneous external magnetic field, besides the inhomogeneous field resulting from 
the magnetic dot, brings a qualitative difference in the free energy diagram. For values of the external magnetic 
field between the average field of the dot in the positive and the negative region, we observe a strong interplay of 
the different superconducting giant states. Moreover, with increase of the thickness of the dot, we obtain two strong 
minima in the free energy, when the total field is approximately equal to zero in the p < Ri, and the p > Ri region, 
respectively, where i?i is the radius of the dot. 

Here, we must emphasize that our results are only valid in the limit of very thin disks. It allowed us to separate 
the two GL equations and thus we neglected the magnetic field created by the superconducting currents. 
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APPENDIX: Eigenfunctions of the linearized first GL equation in the presence of an 

inhomogeneous magnetic field 
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A. Disk magnetic field profile 

The vector potential distribution is determined by the piecewise function 



Hp) 



Hop/2, 

— HoutP/2 + HoutR2l2p, 

0, 



0<p<i?i, 
Ri<P< R2, 
R2< P<R- 
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{III) 
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The eigenfunctions of Eq. (h5f) are expressed in the following way 

fL,n (p) = 



f'ip), (r) 

bif('ip)+b2fi'ip), {II) 

d^fl''{p) + d2fi^'{p), {III) 



(A.2) 
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J„i (a;) and y,„ (a;) are the Bessel functions of the first and second kind, M (a, c, y) and t/ (a, c, j/) are the Kummer 
functions. To find the unknown constants 61(2), '^1(2) and the eigenvalue A we have to join the different parts of 
/i.n (p) and its derivatives at i?i and R2 as well as to use the boundary condition {df/dp)\j^ = 0. 

To simplify the next calculation let us introduce the following notations: Ki — /^(i?i), /i = //^ (i?i), A'^i — 



fi' (i?i), /3 - fC (i?2), A^3 - fi' (^2), P3 = f(" {R2), O3 = fi'' (i?2), A'2 
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have the simple form 



With these notations the system of equations which determine the coefficients &i(2), c^i(2) and A 



From the first four equations of the system (| 



61/1 + b2Ni = 

61/2 + 62iV2 = 

bih + b2N3 - diP3 - d2Q3 - 

bih + b2Ni-diPi-d2Q4-- 

diP5+d2Q5 -- 


= Ki, 

= K2, 

= 0, 
= 0, 
= 0. 


stem (A.3) we obtain: 




^ _ A'iA^2 - i^2iVi 
^ I1N2-I2N1 ' 
I1K2 - I2K1 
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rfl = 



{K1N2 ~ K2N1) (/3Q4 - Q3/4) + {hK2 - I2K1) {Q4N3 - Q3N4) 

P3Q4 - P4Q3 

{K1N2 ~ K2N1) {Psh - hPi) + {hK2 - hKi) {NiPs - N3P4) 



'^^ PsQi - P4Q3 

Inserting these results into the last equation of system (|A.3| ) results into the non-linear equation for A: 
P5 [{K1N2 - K2N1) (/3Q4 - Q3/4) + {hK2 - I2K1) [QiNs - QsN^)] 



(A.5) 



-Q5 [{K1N2 - i^2A^i) (P3/4 - HPa) + {hK2 - hKi) (N^P^ - N^Pi 



0. 



To obtain the correct A values we have to exclude from the spectrum of solutions of Eq. ( A.5 ) those which results in 
zeros of both denominators I1N2 — I2N1 and P^Qa, — PiQz- 

In the large radius limit R — > 00, i.e. which is equivalent to the thin film limit, the Bessel functions in the (Ill)-region 
have to be replaced by their asymptotics. Substituting them in the equation diP^ + d2Q5 = we obtain 



A 



^'^W 



-(3 + 2|L|) -arctan — 

4 0.2 



fL,n (p) near the edge of the sample is equal to i/df+df. 



(A.6) 



B. Ring magnetic field profile 



The vector potential distribution in this case is 

0, 



Hop/2 - HoRy2p, 



HoiitP/2 + HoutR2/2p, 



0, 



0<p<Rd, (I) 

Rd<P<Ri, {II) 

Ri<P<R2, {III) 

i?2 < P < i?, {IV) 



and, therefore, the eigenfunctions of Eq. (IT^) are expressed as follows: 

/i,n {p) = 



where 



with 



f'{p). {I) 

bif('{p) + b2.fi'{p), {II) 

difl''{p) + d2fi''{p), {III) 

ei/r(p)+e2/|^(p), {IV) 



f {p) = fr {P) = J\L\ (VT+Ap) , 

f(' {p) = (ffopV2)'^"'^' exp i-H.pyA) M {-,.,, jj (A) , \L}j\ + 1, HopV2) , 

f" {P) - (ffopV2)"^"'^' exp (-ffopV4) U {-v^ji (A) , \L}j\ + 1, Hop" 12) , 

fi'' {p) = {HoutpV2)^'''"'^'^ exp {-HoutpV^) M {-v^ju (A) , \L*ui\ + 1, i^o.tpVs) 

/f ^ {p) = (ifo„ipV2)l''"^l/' exp (-iJo„,pV4) V {-vnjii (A) , |il,,| + 1, i/o«tpV2) 

H"" {P) = y\L\ (x/T+Ap) , 



(A.7) 



(A., 



T* 
T* 



t^n.II (A) = - 



L + H^Rl/2 
L- 
1 



IIoutR2/2, 



II 



Vn,III (A) = - 



I + IL 



iii\ 



L 



III 



1 + A 

2iJo ' 

1 + A 

2Hnut 
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The unknown constants ^1(2), £^1(2)7 61(2) and the eigenvalue A are found by joining the different parts of /l_„ (p) and 
its derivatives at Rd,Ri and R2 as well as from the boundary condition [df /dp)\ j^ = 0. The corresponding system 
of equations has the form: 

biIi+b2Ni ^Ki, 

bil2+b2N2 ^ K2, 
blh + 62iV3 - diPs - ^2^3 = 0, 

61/4 + b2Ni ~ diPi - d2Qi = 0, (A.9) 

diPs + d2Q5 - eiSz - e2T^ = 0, 
diPQ + d2Qe - ciSq - e2Te = 0, 
eiiSy + e2Tj = 0. 



where the following notations were introducediiiTi = f^ {Rd), Ii — fi^ (Rd), Ni — fj' (Rd), K2 — 

, ^3 = fl' {Ri), N3 = fi' (i?i), P3 = //" (i?i), Q3 = fi'' (Ri), h 



df'ip) 



dp 



dfi'ip) 
dp 

Na = 



, A^2 = 

p=Rd 
9/2" (p) 



dfl'(p) 



dp 



P=Ri 



dp 

Pi 



p=Rd 



p=Rd 



_ 9fi'(p) 
- dp 



P=Ri 



T5 = /i^(i?2), Pe 



T7 



dfi^(p) 
dp 



P=R2 



P=Ri 



Qi = 

dfV'(p) 
dp 



dfi"{p) 



dp 



P=R2 



P=Ri 



P. 



/r(i?2), Q, = fi"{R2), S, = /r(i?2), 



P=R2 



n 



df^iP) 



P=R2 



Sj 



dp 



p=R 



p=R 



The constants &i(2) and ^1(2) are determined by the expressions (A.4) and 



ei 



_ (P^Te - nPe) d, + jQ^Te - T^Qe) ^2 



SsTq — SrTi 



'6^5 
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{S5P6 - PsSe) di + jS^Qe - Qs^e) ^2 
SsTg — SrTc 



6-^5 



From the last equation of system (A.9) 



[{K1N2 - K2N1) {hQi - Q3I4) + {I1K2 - I2K1) {Q4m - Qsm)] 

X [Sr (PsTe - T^Pe) + T7 {S^Pe - P^Se)] 

+ [{KiN2 - K2N1) (P^h ~ hPi) + {hK2 ~ I2K1) (NiPa - N^Pi)] 

X [^7 (Q5T6 - TsQe) + Tr {S^Qe - Q^Se)] = 



(A.IO) 



we obtain the spectrum of A values (notice, that we have to exclude the zeros of the denominators /iA'2 — I2N1, 
P3Q4 - PaQz and S'sTg - SqT^). 

In the large radius limit R —f 00, similar considerations have to be made as for the previous case of the disk profile. 
/l,™ (p) near the sample edge is equal to \/^+e^ and 



A = -l 



1 
i?2 



-(3 + 2|i|) -arctan — 

4 62 



(A.ll) 
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FIG. 1. The magnetic field profile of the magnetic disk (solid curve) and the corresponding model step profile (dashed curve). 

FIG. 2. The magnetic field dependence of the lowest eigenvalues of the linearized first GL equation for different angular 
momenta L and for (a) Ri/£, = 1.5, (b) Ri/^ = 3.0, and (c) Ri/(, = 4.5. The corresponding curves for the first radial excited 
state, i.e. n = 1, are presented in figures (d-f). The top axis gives the flux through the positive magnetic field region in units 
of the flux quantum. 

FIG. 3. The Cooper pair density for the giant vortex states with angular momenta L = 0, 1, 2 (sohd curves), the antivortex 
states L = —1, —2 (dotted curves) and the ring vortex state n = 1, L = (dashed curve) at the magnetic field Hin ~ 0.75Hc2 
for (a) Ri/^ = 1.5, (b) Ri/^ = 3.0, and (c) Ri/^ = 4.5. 

FIG. 4. The free energy of the giant vortex states with different angular momenta L as a function of the applied magnetic 
field in the positive region for (a) Ri/^ — 1.5, (b) Ri/S, ~ 3.0, and (c) Ri/^ = 4.5. Only the states with \L\ < 10 are shown. 
The insets depict the magnetic field dependence of the disk magnetization for the ground giant vortex state for different values 
of the Ri/^ parameter. Figure (d) gives the positive flux captured in the superconducting disk for the different L- ground states 
as the ground state as function of Ri/S,. 

FIG. 5. The Cooper pair density for the giant vortex states with angular momenta L = 0, 1, 2 (solid curves), the anti-vortex 
states L — —1, —2 (dotted curves) and the ring-shaped vortex state n = 1, L = (dashed curve) at the magnetic field 
Hin = 0.75Hc2 for (a) R/^ = 9.0, and (b) R/^ = 12.0. 

FIG. 6. The Cooper pair density for (a) the Meissner state (L = 0), and (b) the ring-shaped vortex state n = 1, L = at 
the magnetic field Hi„ = 0.75Hc2 for different sizes of the superconducting disk. 

FIG. 7. The free energy of the giant vortex states with different angular momenta L as a function of the applied magnetic 
field for (a) R/^ = 9.0, (b) R/^ = 12.0, (c) R/^ = 18.0, and (d) R/^ = 20.0, with Ri/^ = 4.5, and 7?2/? = 6.0. Dashed curves 
represent the energy of the anti- vortex states. 
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FIG. 8. The magnetic field dependence of tfie disk magnetization for the ground giant vortex state corresponding to the 
states in Fig. §, for (a) R/^ = 9.0, (b) R/S, = 12.0, and (c) R/S, = 18.0. 

FIG. 9. The free energy of the giant vortex states with different angular momenta L as a function of the external magnetic 
field for (a) Ri/^ = 5.25, i?2/C = R/(, = 6.0, and (b) Ri/^ = 4.5, i?2/C = 6.0, R/^ = 9.0. Dashed curves illustrate the energy 
of the anti- vortex states. In (a) dotted curves represent the giant-giant multivortex states while in (b) dotted curves depict the 
free energy of the giant-ring multi-vortex states, and the solid curves correspond to the giant-giant states. Thick solid curve 
illustrates the energy of the ring vortex state, i.e. (n, L) — (1, 0) 

FIG. 10. Transition between two giant vortex states shown through the contour plots of the superconducting wave function 
density for the (0 : 4) giant-giant multi-vortex state for different values of the magnetic field in the positive field region 
(corresponding to the magnetization of the dot). 

FIG. 11. Contour plots of the superconducting density for the ground state and the corresponding phase contour plots (see 
Fig. H(b)) for different values of the magnetic field in the positive region. 

FIG. 12. Contourplot of the superconducting density for the ((5,0):(0, 1)) giant-ring multi-vortex state (see Fig. ^(b)) for 
different values of the applied magnetic field. 

FIG. 13. (a-c) Contour plots of the superconducting density for the (0:7), (0:8), and, (0:9) giant-multi anti-vortex states, 
respectively, for different values of the magnetic field Hin; (d-f) the corresponding contour plots of the phase of the supercon- 
ducting wave function density. Notice that the phase near the boundary is near or 27r but due to the finite numerical accuracy 
it oscillates between 2-k — e and 2tt + e where ^ ~ 10~^. 

FIG. 14. The free energy of the giant vortex states with different angular momenta L as a function of the external magnetic 
field for Rd = 2.1^, Ri = 6.6^, R2 ~ 8.1^, R/^ — 9.0. Dashed curves depict the anti-vortex states and dotted curves represent 
the free energy of the giant-giant multi- vortex states. A bird view of the magnetic field profile is given in the inset. 

FIG. 15. The Hin — Ri equilibrium vortex phase diagram for a thin superconducting disk with Rd/(, ~ 0.0, -R2/C ~ 6.0, 
and i?/^ = 6.0. Dashed curves indicate transitions between different giant vortex states and the thick shaded area denotes the 
multi-vortex region. The normal/superconducting state transition is given by the thin solid curve. 

FIG. 16. The Hin — R phase diagram for the ground state of a thin superconducting disk with Rd/^ ~ 0.0, Ri/S, = 4.5, and 
^2/C = 6.0. Solid curves indicate transitions between different vortex states including multi- vortex (giant-multi anti- vortex 
and ring-giant) regions. Dashed curves denote the transitions between different metastable giant vortex states. 

FIG. 17. The Hin — Rd phase diagram for the ground state of a thin superconducting disk with Ri — Rd — 4.5^ and 
R2 — Rd = 6.0^, i.e. a ring inhomogeneous magnetic field distribution. The same curve convention is used as in Fig. [16l 

FIG. 18. The configuration: a superconducting disk with radius R and thickness d with a ferromagnetic dot with radius Ri 
and thickness dd, which is placed on top of it. 

FIG. 19. The free energy of different vortex states as a function of the magnetic moment of the ferromagnetic dot with 
(a) Ri/^ = 4.5, i?2/C = R/^ = 6.0, and (b) Ri/^ = 4.5, i?2/C = R/^ = 9.0. In (a), dotted curves depict the free energy of the 
giant-giant multi-vortex states while in (b) dotted curves denote giant-ring states. The bold curve gives the energy of the ring 
vortex state, i.e. (n, L) — (1,0). The energy of the anti-vortex states is given by the dashed curves. 
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FIG. 20. The magnetic field profile as produced by a current loop (solid curve) and its corresponding model profile (dashed 
curve) . 

FIG. 21. The free energy of different vortex states as a function of the current in a loop with (a) Ri/(, = 7.5, 
R2/(. = R/S, ~ 9.0. Dashed curves depict the free energy of the anti-vortex states and dotted curves denote the giant-giant 
multi- vortex states. 

FIG. 22. The free energy of giant vortex states as a function of the external magnetic field in a superconducting disk with 
magnetic dot on top of it with parameters Ri/(, — 4.0, R/S, = 6.0 and (a) dd = 0.1^, (b) dd — 0.5^, and (c) dd = 1.0^. Dashed 
curves depict the free energy of the anti-vortex states. Insets show the magnetic field profile inside the superconducting disk 
created by the magnetic dot. 

FIG. 23. The free energy of giant vortex states as a function of the external magnetic field in a superconducting disk with 
current loop on top of it with parameters Ri/(, = 7.5, R/(, — 9.0 and (a) I/Io = 5.0, (b) I/Io = 10.0, (c) J/Jo ~ 15.0, and (d) 
//Jo ~ 20.0. Dashed curves depict the free energy of the anti-vortex states. Insets show the magnetic field profile inside the 
superconducting disk created by the current loop. 
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